The Big-Bang Singularity 
in the framework of a Generalized Uncertainty Principle 
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We analyze the quantum dynamics of the Friedmann-Robertson-Walker Universe in the context 
of a Generalized Uncertainty Principle. Since the isotropic Universe dynamics resembles that of a 
one-dimensional particle, we quantize it with the commutation relations associated to an extended 
formulation of the Heisenberg algebra. The evolution of the system is described in terms of a massless 
scalar field taken as a relational time. We construct suitable wave packets and analyze their dynamics 
from a quasi-classical region to the initial singularity. The appearance of a non singular dynamics 
comes out as far as the behavior of the probability density is investigated. Furthermore, reliable 
indications arise about the absence of a Big-Bounce, as predicted in recent issues of Loop Quantum 
Cosmology. 

PACS numbers: 98.80.Qc;11.10.Nx 



I. INTRODUCTION 

The classical theory of gravity, Einstein's General Rel- 
ativity (GR), is very well tested on macroscopic scale, i.e. 
it is expected to work up to nuclear density. On the other 
hand this theory implies the very well-known singularity 
theorems The prediction of singularities undoubt- 
edly represents a breakdown of GR in that its classical 
description cannot be expected to remain valid at the 
extreme condition near a spacetime singularity. One of 
the most important example is the Big-Bang singularity 
appearing in cosmological models. 

In this work we analyze the quantum features of the 
cosmological singularity in the framework of minimal 
length uncertainty relations. The existence of a fun- 
damental minimal scale is long expected in a quantum 
theory of gravity and it invites the possibility that there 
are corrections to the usual Heisenberg uncertainty prin- 
ciple such it becomes what is known as a Generalized 
Uncertainty Principle (GUP). Interest in minimal length 
or Generalized Uncertainty Principle has been motivated 
by study on perturbative Sting Theory [2 L considerations 
regarding the proprieties of black holes 3[ and de Sitter 
space However, in recent years, a big amount of work 
has been made in this active field in a wide variety of di- 
rections (see for example [f| and the references therein) 
and here we present the first application to a minisuper- 
space dynamics. 

In view applying the idea a cutoff scale to the min- 
isuperspace, the GUP approach appears the best physi- 
cally grounded. In fact, it relies on a modification of the 
canonical prescriptions for quantization and, in this re- 
spect, it can be reliable applied to any dynamical system. 
Thus, GUP formalism allow us to implement some pecu- 
liar features of String Theory in quantizing a cosmolog- 
ical model. The application to the cosmological minisu- 
perspace of other noncommutative approaches (like the 



k-Minkowski formulation [6|) would require to address 
its Minkoskian nature as a physical one. But no reliable 
evidence arises that this Lorentian signature has a deep 
physical interpretations. 

The studied model is the isotropic flat Fricdmann- 
Robertson- Walker (FRW) (k = 0) cosmological model 
with a massless scalar field. By the minisuperspace re- 
duction we are able to treat the scale factor as a coor- 
dinate of the particle (Universe), i.e. in the framework 
of GUP a nonzero minimal uncertainty in the isotropic 
volume of the Universe appears. Despite the simplicity of 
the model, this deserves interest for mainly two reasons: 
i) the scalar field can be regarded as an "emergent time" 
for the quantum evolution, ii) in the Wheeler-DeWitt 
framework the unavoidable classical singularity cannot 
be solved and the wave packet follow a classical trajec- 
tory up to the "initial" singularity @, Q • In the GUP 
approach, as we will see, the situation is very different. 
In fact the wave packet, peaked at late times (at energies 
much smaller than the Planck's one), "escape" from the 
classical trajectory in the dynamics toward the cosmolog- 
ical singularity. Therefore the probability density to find 
the Universe near the classic time where the singularity 
appears goes to zero and, in some sense, our quantum 
Universe approach stationary states "near the Planck- 
ian region" . In this sense the cosmological singularity is 
solved by the modified Heisenberg algebra. 

The structure of the letter is the following. In Sec- 
tion II and III the quantum-mechanical structure which 
underlines the GUP framework and the canonical quan- 
tization of the FRW (fc = 0) cosmological model with 
a massless scalar field are briefly reviewed, respectively. 
Section IV is devoted to a discussion on the scalar field 
taken as a relational time. In Section V the "generalized 
quantization" of the model is performed and either the 
probability density to find the Universe near the classical 
singularity, either the wave packets dynamics are investi- 
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gated in detail. In Section VI concluding remarks follow. 



II. QUANTUM MECHANICS WITH GUP 



This is nothing but a generalized Fourier transformation, 
where in the limit f3 — > the ordinary position wave 
function ip{C,) = (CIV') is recovered. 



In this section we briefly review some aspects and 
results of a nonrelativistic quantum mechanics with 
nonzero minimal uncertainties in position [9j. We con- 
sider the Heisenberg algebra generated by q and p obey- 
ing the commutation relation (we work with h = c = 1 
units) 



[q,p] =z(l + /3p 2 



(1) 



where (3 is a "cutoff" constant. This commutation rela- 
tion leads to the uncertainty relation 



AqAp > 1 (1 



- f3(Ap) 2 + p(p) 



(2) 



which appears in perturbative String Theory [2j. This 
Generalized Uncertainty Principle implies an absolutely 
minimal uncertainty in position (considering (p) = 0) 
Ago = \fp- The existence of a nonzero uncertainty in 
position implies that there cannot by any physical state 
which is a position eigenstate since an eigenstate would 
of course have zero uncertainty in position. This consid- 
eration force us to work on the momentum space where 
we have the following representation: 

pip{p) = ptp(p), qtp(p) = i(l + (ip 2 )d p ip{p). (3) 

In order to recover information on positions we have to 
study the states which realize the maximally allowed lo- 
calization. Such states \4>T ) of maximal localization 
around a position ( have the proprieties (V'™ / |q|' ( /'™ / ) = C 



and (AgW 



Ago and obey the following equation 



W) = o, 



(4) 



where we remember that the absolutely maximal local- 
ization can only be obtained for (p) = 0. Thus, the 
maximally localized states in the momentum space read 
explicitly as 



i,f{p) 



1 



(1 +j9p 2)l/2 1 ^ )> 



exp 



(5) 



which, in the limit (3 — > 0, reduce to the ordinary plane 
waves. At this point we can project an arbitrary state \ip) 
on the maximally localized states \tp 1 Q l ), in order to ob- 
tain the probability amplitude for a particle being maxi- 
mally localized around the position £ (i.e. with standard 
deviation Ago). We call these projections the "quasipo- 
sition wave function" -0(0 = (V ; 7 l/ |V'); explicitly we have 



dp 



(l + /3p 2 )3/ 2 



ex P I * 7^ ) W)- 



VP 



(6) 



III. CANONICAL QUANTIZATION OF THE 
FRW MODEL 

The canonical quantization of the homogeneous, 
isotropic, flat (k — 0) cosmological model with a massless 
scalar field is reviewed (for more details see The 
Hamiltonian constraint for this model has the form 



-9np 2 X 



8tt x 



x = a 3 , 



(7) 



8nlp is the Einstein constant and 



where k = 8itG 
a is the scale factor. In the classical theory, the phase 
space is 4-dimensional, with coordinates (x,p x ; 4>,p<p). At 
x = the physical volume of the Universe goes to zero 
and the singularity appears. Since 4> does not enter the 
expression of the constraint, is a constant of motion 
and therefore each classical trajectory can be specified in 
the (cc, 0)-plane. Thus <f> can be considered as an internal 
time (see below) and the dynamical trajectory reads as 



±- 



1 



"V2l? 



(8) 



where xq and 4>$ are integration constants. In this equa- 
tion, the plus sign describes an expanding Universe from 
the Big-Bang, while the minus sign a contracting one 
into the Big-Crunch. We now stress that the classical 
cosmological singularity is reached at <f> = ±00 and every 
classical solution, in this model, reaches the singularity. 

At quantum level the Wheeler-DeWitt equation, asso- 
ciated to the constraint ([7|), tells us how the wave func- 
tion ^{x, (f>) evolves as (f> changes; in this respect we can 
regard the argument <fi of ^(x, 4>) as an "emergent time" 
and the scale factor as the real physical variable. In or- 
der to have an explicit Hilbert space, we perform the 
natural decomposition of the solution into positive and 
negative frequency parts. Therefore, the solution of this 
Wheeler-DeWitt equation has the very well-known form 



-1/2 



(Aa 



-27 



Bx l ^\ 



(9) 



where 7 = ^iAf 2 — l) 1 / 2 > and e 2 being the eigen- 
value of the operator S/247TK defined below. Thus the 
spectrum is purely continuous and covers the interval 
(•\/3/2Zp, 00) 0]. The wave function ^(x, <fi) is of pos- 
itive frequency with respect to the internal time 4> and 
satisfies the positive frequency (square root) of the quan- 
tum constraint ([7]); we deal with a Schodinger-like equa- 
tion —id,),^ = v / S v I / , where S = 2Air Kxp" x 2 x . 

In order to examine the behavior of the classical sin- 
gularity at quantum level we have to clarify a general 
criteria for determining whether the quantized models 
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actually collapse [10(. Unfortunately there is not such 
a rigorous criteria yet. An early idea was to impose the 
condition that the wave function vanishes at the singular- 
ity a = ll|, but this boundary conditions has little to 



do with the quantum singularity avoidance. It seems bet- 
ter to study the expectation values of observables which 
classically vanish at the singularity. In fact, Y$>(a = 0, t)\ 2 
is merely a probability density and thus, for example, one 
might have an evolving state that "bounces" (i.e. a non- 
singular wave packet), even if |\&(a = 0, t)\ ^ for all 
t. On the other hand, if one could find a wave packet so 
that the probability Pg = \^(a, t)\ 2 da ~ for 5 being 
a very small quantity, then one could reasonably claim 
to have a nocollapse situation. 

Let's now come back to the canonical FRW model. It 
is not difficult to see that, in this framework, the unavoid- 
able classical singularity is not tamed by quantum effects. 
In fact, if one starts with a state localized at some initial 
time, then its peak moves along the classical trajectory 
and falls into the classical singularity. Additionally, from 
the eigenfunctions ((9]) it is clear that the probability de- 
fined above diverges indicating that the Wheeler-DeWitt 
formalism does not solve the classical singularity. 



IV. ON THE SCALAR FIELD AS A 
RELATIONAL TIME 



Let us now discuss in some details the role of a mat- 
ter field, in particular of a scalar field (f>, as definition 
of time for the quantum dynamics of the gravitational 
field. As well-known, in any diffeomorphism-invariant 
quantum field theory (like Quantum General Relativity) , 
the Schodinger equation is replaced by a Wheeler-DeWitt 
equation, in which the time coordinate disappears from 
the formalism. This is the so-called problem of time and 
it is one of the major conceptual puzzle in Quantum 
Gravity (for an exhaustive discussion on it see [13]). To 
overcome this difficulty it was argued that the "physical 
time" has to be described using some of the canonical 
variables of the gravitational field. In particular it is 
possible to distinguish between an intrinsic time, which 
is constructed entirely from the 3-metric, and an extrin- 
sic time, which is constructed from both the metric and 
the extrinsic curvature. However, it was shown how it is 
difficult to provide a satisfactory definition of time using 
both these concepts. Moreover, a physical clock has to be 
made of matter with certain definite proprieties. There- 
fore it seems natural to use a matter field to define a 
relational time. The concept of a relational time is prin- 
cipally based on the idea about the absence of the time 
at a fundamental level. Thus, the dynamics of a system 
can be described only on a relational point of view, i.e. 
a field has to be evolved with respect to another one (for 
more details on this approach see [HI]). From this per- 
spective, we will evolve the physical degrees of freedom 



of the gravitational field toward a massless scalar field 
and in this respect it will be treated as an external time. 

In particular, in quantum cosmology, the choice of a 
scalar field as an external time appears as the most nat- 
ural one. In fact, a monotonic behavior near the classical 
singularity of 4> as function of the scale factor (more pre- 
cisely the variable which describe the isotropic expansion 
of the Universe) always appears. Let us consider the case 
of the Bianchi IX model, i.e. the most general (together 
with Bianchi VIII) homogeneous cosmological model, in 
the presence of a massless scalar field <j>. In such a case, 
the Hamiltonian constraint Hgray + H^ w 0, has the form 



4k 



8tt a 3 



0, (10) 



where the a variable describes the isotropic Universe ex- 
pansion, f3± are associated to the space anisotropics and 
the potential term V((3±) accounts for the spatial curva- 
ture of the model. Since the determinant of the 3-metric 
is given by h = a 3 , it is easy to recognize that the clas- 
sical singularity appears for a — ► 0. Therefore near the 
cosmological singularity i.e. in a purely quantum region, 
the potential term can be neglected (for a discussion on 
the consistency condition which ensures that the quasi- 
classical limit of the Universe dynamics is reached before 
the potential term becomes important see [3|)- As be- 
fore, we will regard the scalar field as the time for the 
dynamics and, in this respect, the equations of motion of 
the system read as 



a(<p) = B exp 



A<j> 



VA 2 



, Va{4>) 



A 
B 



exp 



A4> 



VA 2 



(11) 

where A and B are integration constants and c = p\ + 
p_ = const (p + — const and p_ = const) . 

An important point has now to be stressed. As well- 
known [15| , the dynamics toward the cosmological singu- 
larity of a generic inhomogeneous Universe is described 
point by point, by that one of a Bianchi IX model. More 
precisely, it is possible to show that the spatial point, 
toward the singularity dynamically decouple and there- 
fore the spatial geometry can be viewed as a collection 
of small patches, each of which evolves independently as 
a Bianchi IX model. In this respect, the above mono- 
tonic relation between a massless scalar field and the 
isotropic variable of the Universe, is a real general fea- 
ture of the gravitational field and this clarify the choice 
of a scalar field as a relational time. Therefore we can 
conclude that, not only in the symmetric model treated 
in this work but also in the most general cosmological 
model, the scalar field can be regarded as a good time 
for the gravitational dynamics. 
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V. GENERALIZED QUANTIZATION OF THE 
FRW MODEL 



We will analyze the quantization of the FRW (k = 0) 
model in the framework of minimal length uncertainty 
relation. As in the canonical case, let us decompose the 
solution of the "generalized Wheeler-DeWitt equation" 
into positive and negative frequency parts and focus on 
the positive frequency sector. Remembering that we have 
to work in the momentum space, the wave function reads 
as: *ff(p, (f>) — ^E'(p)e 4%/24Tr ' ce *, where from now on p = p x . 
Thus, as soon as we regard the scalar field as an "emer- 
gent time" for the quantum evolution, then we treat in 
the "generalized" way only the real degree of freedom 
of the problem: the isotropic volume x. Therefore the 
quantum equation relative to the Hamiltonian constraint 
(O, considering the representation ([3]), is the following 

M 2 (l+M 2 ) 2 — + 2 A i(l+/i 2 )(l + 2/r 2 ) — +e 2 * = 0, (12) 



now construct states peaked at late times 



df i 



where we have defined a dimensionless parameter /i = 
\fpp. In order to integrate the above equation we intro- 
duce the variable p = tan -1 n, which maps the region 
0</i<ooto0<p< 7r/2. Then, performing an- 
other change of variables: £ = ln(sinp) (— oo < £ < 0), 
equation (fi"2"|) reduces to 

^ + 2 (l^M +eV ^°' (13) 

which can be explicitly solved and whose general solution 
reads as 



* e (0 = c ie -«( 1+Q > (i 



b+e^+^e^ 1 -^ (l + 6_e 2 «). 

(14) 

where a = \/l — e 2 and b± — 1 ± a/(l =F a). At this 
point we have to analyze the "quasiposition wave func- 
tion" relative to this problem in order to make a first 
comparison with the canonical case, in particular with 
the wave function ([9]) . In agreement with formula ([6]) we 
have 



*e(C) = / d£ cxp U + < tan" 

J —OO \ 



Cie -«(i+«) (i + 6+e 2?) + c 2 e-^-^ (1 + 6_e 2 «) 

(15) 

where £, in this case, is expressed in units of V/3- Thus we 
can easily see that our "quasiposition wave function" , i.e. 
the probability amplitude for the particle (Universe) be- 
ing maximally localized around the position £, is nondi- 
verging for all £, as soon as we take the condition C\ = 0. 
We stress that the canonical wave function (the function 
dH)) is diverging at the classical singularity x = 0. 

To get a better feeling with our quantum Universe we 
construct and examine the motion of wave packets. Let's 



*(C,*) 



de fl (e)* e (0e* 



(16) 



where we have defined the dimensionless time t — 
V24 TTKip. In the following we take 17(e) to be a Gaussian 
distribution peaked at some e* <C 1, which corresponds 
to be peaked at energy much less then the Plank energy 
l/lp (we recall that e ~ 0(e lp), where e have dimension 
1 in energy). The analytic computation of the integral 
(|16|) for the wave function (|15p is impossible to perform. 
So, in order to describe the motion of wave packets we 
have to evaluate (fl6|) numerically. 

At first, we want to analyze the most interesting re- 
gion, i.e. where C — 0, which corresponds to the purely 
quantum region, where the physical volume is Planck- 
ian. In fact, if we put (3 ~ 0(l P ), the minimal uncer- 
tainty in position is of order of the Planckian volume. 
The "quasiposition wave function" (|15[) can be expanded 
in order to give the probability density around ( ~ 0: 
|*(C,i)| 2 ~ \A(t)\ 2 + ( 2 \B(t)\ 2 . Therefore, starting with 
a state peaked at some e* <C 1, the probability density 
of finding the Universe "around the Planckian region" is 
|A(t)| 2 , where A(t) reads 



W Jo VT^ (3 - e 2 + 3VT^) 



(17) 



We evaluate the above integral numerically for e* = 1CU 3 , 
cr 2 = 1/20 and we take the constant 2C2 = 1. The 
probability density |>l(t)| 2 is very well approximated by 
a Lorentzian function (see Fig. 1). As we can see form 
the picture, this curve is peaked around t = 0. This value 
corresponds to the classical time for which x{t) — xq (in 
(fTT)) we consider to — 0). Thus, for x ~ 0(lp), the prob- 
ability density to find the Universe in a Planckian vol- 
ume is peaked around the corresponding classical time. 
As a matter of fact this probability density vanishes for 
t — > —00, where the classical singularity appears. This is 
the meaning when we claim that the classical cosmolog- 
ical singularity is solved by this model. 

In order to describe the motion of the wave packet 
we evaluate ^(d t)\ 2 from the integral (flTj)) of the wave 
function (|15p . As before, we consider a wave packet ini- 
tially peaked at late times and let it evolve numerically 
"backward in time" . We use the same parameters for 
the integration performed above. The result of the inte- 
gration is that the probability density, at different fixed 
values of £, is very well approximated by a Lorentzian 
function yet. The width of this function remains, actu- 
ally, the same as the states evolves from large C (10 3 ) to 
C = 0. For all fixed t the probability density is well- fixed 
by a Lorentzian function and the width of these function, 
also in this case, remain almost the same during the evo- 
lution. These states are sharply peaked for £ ~ 0(1) 
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FIG. 1: The points represent the result of the numerical in- 
tegration and are fitted by a Lorentzian Lit) = 0.692/(t 2 + 
31.564) heaving width, at the inflection point, 3.243. 



(which in our units correspond to C ~ 0(yf]3) ~ 0(lp)). 
The peaks of Lorentzian functions, at different £ values, 
move along the classically expanding trajectory ([8j) for 
values of £ larger then ~ 4. Near the Planckian region, 
i.e. when £ S [0, 4], we observe a modification of the tra- 
jectory of the peaks. In fact they follow a power-law up 
to £ = 0, reached in a finite time interval and "escape" 
from the classical trajectory toward the classical singu- 
larity (see Fig. 2). The peaks of the Lorentzian at fixed 
time t, evolves very slowly remaining close to the Planck- 
ian region. Such behavior outlines that the Universe has 
a stationary approach to the cutoff volume, accordingly 
to the behavior in Fig. 2. 

This peculiar behavior of our quantum Universe is dif- 
ferent from other approaches to the same problem. In 
fact, recently, it was shown how the classical Big-Bang 
is replaced by a Big-Bounce in the framework of Loop 
Quantum Cosmology (LQC) [§]. Intuitively, one can ex- 
pect that the bounce and so the consequently repulsive 
features of the gravitational field in the Planck regime 
are consequences of a Planckian cut-off length. But this 
is not the case. As matter of fact we can observe from 
Fig. 2 that there is not a bounce for our quantum Uni- 
verse. The main differences between the two approaches 
resides in the quantum modification of the classical tra- 
jectory. In fact, in the LQC framework we observe a 
"quantum bridge" between the expanding and contract- 
ing Universes; in our approach, contrarily, the probabil- 
ity density of finding the Universe reaches the Planckian 
region in a stationary way. 

VI. CONCLUDING REMARKS 

In this work we have shown the effects of a modi- 
fied Heisenberg algebra, which reproduces a GUP as ap- 
peared in studies on String Theory on the cosmolog- 
ical singularity. Two main remarks on our analysis must 



t 




FIG. 2: The peaks of the probability density |*(£,i)j 2 are 
plotted as functions of t and ln(£). The points (resulting from 
numerical computation) are fitted by a logarithm 0.050 ln(£) + 
0.225 for C, > 4 and by a power law 0.067£ 1060 for £ £ [0,4]. 



be stressed: i) we have restricted our investigation to an 
isotropic, homogeneous model, although a quantum Uni- 
verse has to be described by a generic cosmological model 
where all the symmetries are removed; ii) the width of 
the Lorentzian functions in t, and so the probability den- 
sity, are too wide to formulate a precise description of the 
wave packets toward the Planck era. These shortcomings 
of our analysis will be addressed in a more detailed dis- 
cussion on this topics, which will appear elsewhere. 

In this model, the resolution of the classical singularity 
relies on the mathematical framework used. In fact, the 
features of the "quasiposition wave function" are respon- 
sible for the possibility to expand the probability density 
|4 , (£, t)\ 2 near £ ~ 0, while the corresponding probability 
density near a ~ diverges in the Wheeler-DeWitt case. 

Certainly, the main result of the present work is that 
in the framework of GUP no evidence for a Big-Bounce 
seems to come out. In fact the behavior of the wave 
packet near the Planckian region does not allow for a 
bounce, in the sense that the peak of the probability 
density approaches the region £ ~ by a power-law and 
it does not rejoin itself to the contracting region. The 
big width of the probability density is not qualitatively 
a real problem for the singularity absence. In fact, near 
£ ~ 0, the behavior we have found is not the classical 
logarithmic one even taking the error into account. 



ACKNOWLEDGMENTS 

We would like to thank Roberto Guida for his help in 
generating the manuscript plots. 



* Electronic address: battisti@icra.it 



6 



' Electronic address: montani@icra.it 

[1] S. W. Hawking and G. F.R.Ellis, The Large Scale Structure 
of Space-Time (CUP, Cambridge, 1973). 

[2] D.J. Gross and P.F.Mendle, Nucl.Phys.B 303 (1988) 407; 
D.J. Gross, Phys.Rev.Lett. 60 (1988); 1229 D.Amati, 
M.Ciafaloni and G.Veneziano, Phys.Lett.B 216 (1989) 
41; E.Witten, Phys.Today 49 (1997) 24; G.Amelino- 
Camelia et al., Mod.Phys.Lett.A 12 (1997) 2029. 

[3] M.Maggiore, Phys.Lett.B 304 (1993) 65. 

[4] H.S.Snyder, Phys.Rev. 71 (1947) 38. 

[5] A.Kempf and L.Lorenz, Phys.Rev.D 74 (2006) 103517; 
A.Ashoorioon, A.Kempf and R.B.Mann, Phys.Rev.D 
71 (2005) 023503; L.N.Chang et al., Phys.Rev.D 65 
(2002) 125027; F.Brau and F.Buisseret, Phys.Rev.D 74 
(2006) 036002; J.Y.Bang and M.S.Berger, Phys.Rev.D 
74 (2006) 125012. 

[6] S.Majid and H.Ruegg, Phys.Lett.B 334 (1994) 348. 

[7] W.F.Blyth and C.J.Isham, Phys.Rev.D 11 (1975) 768. 



[8] A.Ashtekar, T.Pawlowski and P.Singh, Phys.Rev.Lett. 
96 (2006) 141301; A.Ashtekar, T.Pawlowski and P.Singh, 
Phys.Rev.D 73 (2006) 124038. 
[9] A.Kempf, G.Mangano and R.B.Mann, Phys.Rev.D 52 
(1995) 1108; A.Kempf, J.Math.Phys. 38 (1997) 1347. 

[10] M.J.Gotay and J.Demaret, Phys.Rev.D 28 (1983) 2402. 

[11] B.DeWitt, Phys.Rev. 160 (1967) 1113. 

[12] C.J.Isham, Canonical Quantum Gravity and the Problem 
of Time, arXiv:gr-qc/9202033 K.Kuchar, Time and In- 
terpretation of Quantum Gravity, in Pre 4th Canadian 
Conf, World Scietific Singapore 1992. 

[13] C.Rovelli, Quantum Gravity (CUP, Cambridge, 2004). 

[14] M.V.Battisti and G.Montani, Phys.Lett.B 637 (2006) 
203. 

[15] V.A.Belinski, I.M.Khalatnikov and E.M.Lifshitz, 
Adv.Phys. 31 (1982) 639. 



